Spiral Motion in a Noisy Complex Ginzburg-Landau Equation 
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The response of spiral waves to external perturbations in a stable regime of the two-dimensional 
complex Ginzburg-Landau equation (CGLE) is investigated. It is shown that the spiral core has a 
finite mobility and performs Brownian motion when driven by white noise. Combined with simula- 
tion results, this suggests that defect-free and quasi-frozen states in the noiseless CGLE are unstable 
against free vortex excitation at any non-zero noise strength. 
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Vortices play a central role in the ordering of two- 
dimensional equilibrium systems characterized by a com- 
plex scalar order parameter, e.g., planar magnets and 
Bose-condensates . They also arise in various nonequi- 
librium situations, e.g., the Belouzov-Zhabotinsky reac- 
tion, certain regimes of fluid flows, and the contraction 
of heart muscles |^-||] . Their long-time motion is deter- 
mined by their coupling to the slow modes of the system, 
and hence is expected to acquire a universal character. 

A distinctive but common feature exhibited by a vortex 
in a nonequilibrium oscillating state is the emission of 
spiral waves which change the oscillating frequency of the 
entire system. This behavior is seen, e.g., from solutions 
of the generic model for these situations, the complex 
Ginzburg-Landau equation (CGLE) l2ir 



dta 



(1 + ic)\a\^a + (1 + ib)\7^a, 



(1) 



where a(x, t) is a complex scalar field, and b and c are 
real numbers. Spiral waves are observed for any b ^ c. 
In certain regimes of the parameter space (5, c) , spiral 
defects are spontaneously generated and undergo violent 
motion, while in other cases they are quite static and lock 
into a quasi- frozen structure |6| . 

In this Letter, we derive mobility relations for a single 
spiral in the regime of Eq. (|l|) , where the uniformly oscil- 
lating state as well as spiral wave excitations are linearly 
stable. We show that, in contrast to the equilibrium XY 
model (which corresponds to 6 = c = 0), a vortex defect 
possesses a nonzero mobility and responds only to per- 
turbations which fall within a distance ^, the "screening 
length", from the spiral core. When a noise source is 
present, the above property yields a diffusion constant D 
for the spiral core proportional to the noise strength. The 
diffusion constant D is calculated through a numerical 
scheme with the result well in accord with direct simula- 
tions of a noisy CGLE. Secondly, we present simulation 
results to show that both the defect-free state and the 
frozen-defect state of the noiseless CGLE disappear when 
the (thermal) noise is introduced. These findings are ra- 
tionalized in the context of noise-driven vortex diffusion 



and thermally-activated nucleation of vortex-antivortex 
pairs. We believe the analysis will help toward under- 
standing various other regimes exhibited by the noise- 
less CGLE, in particular, the melting of the frozen-defect 
state into the defect-turbulence state. 

A single-spiral solution to Eq. (|l|) is well documented 
0,11 . For a single-charged vortex centered at xq, the 
solution takes the form, 



,(x, t) = F{r) exp[i{e + ij{r) - cut)], 



(2) 



where r = |x — Xo| and 9 is the polar angle measured 
from the vortex core. Far away from the core, the so- 
lution approaches a plane wave with ip{r) ~ fcr, where 
the asymptotic wavenumber k is related to the rotation 
frequency as w = c+{b — c)k'^. The dependence of fc on & 
and c is known analytically for \b — c\ <^ 1 and \b — c\ ^ 1, 
e.g. k ~ -c~^exp(-7r/|2c|) for 6 = and |c| < 1 0. 

In the following discussion, we shall limit ourselves to 
the case b = and c 7^ 0. This choice is partly mo- 
tivated by a desire for simplicity, but it also offers an 
interesting framework: all plane waves of wavenumber 
k < kc{c) = (3 -I- 2c^)~^/^ are linearly stable in this case, 
including the homogeneous state a = exp{—ict) at fc = 0. 
For c < Ci ~ 1.08, an isolated spiral is stable and has a 
wavevector fc < fee In the case b = c = 0, Eq. (|^), 
in the presence of noise, reduces to the Ginzburg-Landau 
equation for a superfluid (or a planar magnet), allowing 
the crossover between the equilibrium XY model and the 
nonequilibrium CGLE to be explored. In addition, in- 
creasing c, solutions of (|l|) show a transition from the 
quasi-frozen state to the "defect turbulence" state [D, 
which can then be probed in the presence of noise. 

We now consider the response of the spiral to a weak 
additive broad-band noise r/(x, t). 



dta = a - (1 + ic)\a\'^a + \7^a + ?7(x, t). 



(3) 



To the linear order in 77, we can write the perturbed so- 
lution in the form, 



a(x, t) = as(x, + J dx'dt'gs{x, t; x', t')r]{x.' ,t'), 



(4) 



1 



where 5s (x, t; x', t') is the Green's function for the CGLE 
hnearized around as(x, t). Although the exact form of 
gs can only be obtained numerically, it is instructive to 
examine its behavior when x' is far away from the spiral 
core. This can be done as follows. 

In the far field, the spiral wave is very close to a 
planewave with a wavevector k — ke^- Correspondingly, 
the inhomogenous Green's function at large distances 
can be approximated by the Green's function of the 
planewave. The planewave Green's function can be de- 
composed into contributions from two types of modes, 
the amplitude modes with a time constant of order one, 
and the soft phase modes. Only the latter part survives 
at long times, yielding the asymptotic expression, 

exp[i(k • X — Lot)] 

5p(x + xo,i + to;xo,io)'^ 



1/2 
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where and x±_ are parallel and perpendicular compo- 
nents of x with respect to k, and a II — 1 — 2(l+c^)/^^/(l — 
fc^). (Convective instability sets in at Q!|| — 0.) The group 
velocity of the plane- wave is given by Vg ~ Vk^^' = — 2ck. 

Equation implies that a localized perturbation 
spreads diffusively but, at the same time, its center trav- 
els at the group velocity Vg. Equating Vgt with the dif- 
fusion length t^/^, we obtain a length ^ ~ u"^ ~ |cfc|~^. 
This is the decay length of disturbances in the "up- 
stream" direction (i.e., against Vg), beyond which the 
influence becomes exponentially small at any time. Thus 
only perturbations within the distance ^ from the spiral 
core can have a significant influence on its motion. This 
is in contrast with the equilibrium XY model (c = and 
Vg = 0), where the effect of disturbances far away decays 
only algebraically with the distance. 

A quantitative analysis of the spiral motion under the 
perturbation i] can be done by decomposing the correc- 
tion term in (^) into two parts: (i) shape change of the 
spiral wave, and (ii) translation Xq — > Xg and shift of the 
overall phase V'o = V'(O) ~* V'o- Since the unperturbed 
spiral wave is stable, the flrst type of response decays in 
time when the perturbation is switched off. On the other 
hand, the second type of response corresponds to the zero 
modes of the unperturbed spiral solution and hence does 
not decay away. Under an inflnitesimal translation of the 
core position xq = {xq, yo), the spiral solution changes to 
a'g = as +xoWx +yaWy, where Wx and Wy are the transla- 
tional zero modes in the x and y direction, respectively. 
They are solutions of Eq. (|^) linearized around a steady 
spiral, Eq. (|^). In the presence of the noise 77, the equa- 
tions for the translation zero modes assume the form: 



io = Qx, yo = Cy, 



(6) 



where Cx and Cg are projections of 77 on the zero modes 
Wx and Wy, respectively. 



While the form of the functions Wx and Wy can be ob- 
tained directly by differentiating (^), the calculation of 
Ca; and (y is quite difficult. This is because the operator 
C which deffnes the linearized equation around a spiral 
solution is non-hermitian except at c = 0, so the eigen- 
functions Ua, a = x,y, of the adjoint operator £''' are in 
general different from the eigenfunctions Wa- For c 1, 
we have developed a numerical scheme to determine 
Ua, details of which will be reported elsewhere Quite 
generally, Ua{r, 0) decays exponentially at large distances 
from the spiral core [p|. From the projection formula. 



/ (Pxul^i] 

J CPxUaWa 



(7) 



we see that only the part of the perturbation rj sufficiently 
close to the core contributes to Cq, in accord with the 
qualitative analysis developed above. The exponential 
decay of Ua also ensures the convergence of the integral 
on the denominator, yielding a finite "mobility" to the 
spiral core. This contrasts with the c — case, where 
the integral diverges, and the vortex mobility decreases 
logarithmically with the relevant length scale. 

We now apply the above discussion to the special case 
where the external perturbation takes the form of a weak, 
uncorrelated white noise with zero mean and correlators, 

(?7^(x,t)77^(x',t')) -2T5/3^,5(t-t')<5(x-x'), (8) 

where /3, 7 specify real and imaginary parts of rj and T 
characterises the noise strength. Solving Eq. (^, we 
obtain a diffusion law for the vortex core at long times. 



([xo(t)-xo(0)]2} ^4i?t. 



(9) 



For small T, the diffusion constant is D = jiT, with the 
mobility given by 



/ (Pxulux 
1/ d^xuiwxl' 



(10) 



We have evaluated Eq. (|T^) for several different values 
of c using the numerically determined eigenfunctions 
and Wa, and have also performed Langevin simulations 
of the noisy CGLE (||). In the simulations, the core co- 
ordinates of an initially prepared spiral was followed in 
time at various noise strengths T. Brownian motion of 
the core was observed. In Fig. 1(a), we show the sim- 
ulation data for the mean-square deviation of the core 
coordinates versus time for c = 1.4 and various values 
of temperature T. From the slope of the curves, we de- 
termine the diffusion constant D, which is found to vary 
linearly with T at small T. In Fig. 1(b), we plot the 
mobility against the parameter c, where filled squares 
correspond to data obtained from simulations, and open 
circles data from direct calculations using Eq. The 
agreement between the two sets is apparent. 
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1. Noise-driven spiral diffusion: (a) mean-square de- 





FIG. 

placement (r^) = ([xo(t) — xo(0)]^) vs. time for c — 1.4, 
T = 0.001, 0.002, 0.003, and 0.004 (from bottom to top). 
Data obtained from the trajectory, over At ~ 10^, of a spiral 
initially centered in a disk of radius 128 with no-flux boundary 
conditions, (b) mobility fi = D/T vs. c. 

In the special case c — 0, the problem reduces to that 
of vortex diffusion in an equilibrium XY model, which 
has been studied extensively in the past. In this case, 
the linear operator C becomes hermitian, and Ua = Wa- 
For |c| <C 1, we have analyzed the eigenvalue equation for 
Ua in detail, and found that Ua — Wa at distances up to 
r ~ ^ before it switches over to exponential decay at r > 
^. Using this property and the exponential divergence of 
^ for c ^ 0, we obtain 



/i ~ 2c/7r^ 



(11) 



We have checked that this small c behavior is consistent 
with the trend exhibited by the numerical data shown in 
Fig. 0, though the next order term is significant when c 
reaches a value of order one. 

The discussion presented above on the diffusivity of a 
single vortex serves as a useful basis for understanding 
the behavior of the noisy CGLE (^) under more general 
initial conditions. Specifically, we have examined numeri- 
cally the behavior of the system under two representative 
initial conditions: (a) a uniform, defect-free state with 
a = 1, and (b) a random initial state with \a\ <^ 1. (In 
the latter case the system initially develops into a con- 
figuration with a high density of vortices.) It has been 
shown that, in the absence of noise, case (a) is stable and 
the system remains synchronized (i.e., with bound or no 
phase fluctuations) at all times, while case (b) evolves 
into a "quasifrozen-defect" state for c < c* ~ 1.68 and a 
"defect-turbulence" state for c > c* Our simulations 
show that, when noise is present, the distinction disap- 
pears and the system evolves into a single steady-state 
with a finite density of vortices. 

Figure 2 shows the breakdown, under weak noise (T = 
0.002), of an initially-blocked configuration [Fig.||(a)] ob- 
tained at c = 1.4 by simulating the noiseless equation 
with a random initial condition. Clear spots with a black 
center indicate the depression of \a\ at the vortex cores. 
Switching on the noise, we see the diffusive motion of in- 
dividual vortices, creation of new vortices — in this case, 
through noise amplification by a convective instability 



[clear belt in the lower left corner in Fig. 2(b)] — , and 
annihilation of vortex-antivortex pairs. For weaker noise, 
the whole process occurs on larger timescales, reaching 
an asymptotic state with larger spirals. 




(c) (d) 

FIG. 2. Snapshots of \a\ during the breakdown of a blocked 
configuration at c — 1.4. System of linear size L — 128 with 
periodic boundary conditions, (a): blocked state (T — 0). 
(b,c): initial stages (T = 0.002, t ~ 500,700); note the \A\ 
modulations emitted by the cores of the large spirals, (d): in 
the asymptotic state (T — 0.002, t ~ 5000), spirals diffuse; 
their maximal size is fixed by the strength of the convective 
instability and the noise level. 

For c < Ci, the convective instability is absent and 
vortex-antivortex pairs are created solely through a 
"thermal" activation process. (Although this mechanism 
is also present for c > c^, it does not play a dominant 
role when the noise is weak.) Since the amplitude of |a| 
has to be zero at the vortex core, nucleation of a pair 
requires a rare fluctuation in the noise amplitude, unless 
it happens next to the core of an existing vortex. Indeed, 
most of the newly created "thermal" pairs are found next 
to existing vortices in our simulations. When the noise is 
weak, the probability for the creation of a thermal pair 
can be argued to be proportional to an exponential fac- 
tor, exp{—E/T), where the "activation energy" E is ex- 
pected to be lower when next to an existing vortex than 
in other regions of the system. This is borne out by our 
simulation data, which show the exponential dependence 
on T of the ensemble-averaged time r necessary for the 
nucleation of a thermal pair from defect-free, tilted ini- 
tial conditions (Fig. |3^). We also find an almost linear 
variation of E with phase gradient (Fig. ||b). Moreover, 
E is found not to be sensitive to the value of c fill. 
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FIG. 3. (a) Logarithm of the average nucleation time r 
of a thermal pair of vortices vs. 1/T for fiat, tilted initial 
conditions with mean phase gradient k. From top to bottom: 
increasing k, from fc = to fc ~ 0.3. Each point corresponds 
to 500 runs of a system of size L = 128 at c = 1.25. (Even 
though c > Ci, the convective instability has no influence on 
these results.) (b) Nucleation energy E, as measured from 
the small T behavior of (a) (r (x exp(— i?/r)), vs. k. 

In the case of the equiUbrium XY model, below the 
Kosterlitz-Thouless transition temperature Q|, vortices 
of opposite charge attract one another through a loga- 
rithmic Coulomb potential and remain bounded in pairs. 
For c 7^ 0, however, this attraction is cut off at the dis- 
tance ^, beyond which the interaction becomes exponen- 
tially weak Thus, vortices become "free" when their 
separation is larger than ^. 

For c ^ the diffusive motion of individual vortices 
makes it possible for the nucleated vortices to break away 
from one another. (Most of them are short-lived, but a 
fraction do become "free" On the other hand, it 

also enables two distant, oppositely-charged vortices to 
move into the interaction range ^ and annihilate each 
other. The steady-state density n of free vortices is thus 
controlled by both the "pumping" rate p of creating vor- 
tex pairs of size ^, and the diffusion constant D which 
limits the rate of annihilation. On a mean-field level, we 
can write: p = DrP' . When n is small, we may write 
p{n) = po + Pin^ where po is the nucleation rate of pairs 
in the far field, while pi is the nucleation rate next to a 
free vortex. The resulting equation has a unique positive 
solution, corresponding to a unique steady-state. Our 
numerical experiments agree with this conclusion: from 
both synchronized and random initial conditions, a single 
steady-state seems to be reached JTT[ . 

For c > Ci, the asymptotic state possesses a nonuni- 
form density of vortices. There are two relevant length 
scales in this case: the "noise" length ^ log(l/T) given 
by noise amplification condition, which determines the 
size of the large spirals, and the length scale ^ below 
which vortices of opposite charge annihilate each other. 
(At such large values of c, ^ becomes comparable to the 
size of the vortex core.) The system is better described 
as a gas of large, well-developed spirals with a texture of 
precipitated vortices at their boundaries [Fig. §(d)]. 

In conclusion, we have shown that vortices in the 
CGLE possess a finite mobility and diffuse under the in- 



fluence of an external white noise. The diffusion constant 
at small noise strength is calculated through a numerical 
scheme which compares favorably with direct simulations 
of a noisy CGLE. Our results indicate that the previ- 
ously observed frozen-defect state of the noiseless CGLE 
is unstable against noise, and hence it does not possess 
the kind of metastability exhibited by a typical glassy 
state. Noise also leads to the creation of free vortices 
through either an activated process or a dynamic insta- 
bility. The steady-state of the system becomes unique 
and contains a finite density of free vortices, in contrast 
to the behavior of the equilibrium XY model below the 
Kosterlitz-Thouless temperature. 
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